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The neutralino pair annihilation into two photons in our galactic halo gives a 
robust dark matter signal, since it would give a quasi-monotonic gamma ray. This 

•rH , 

^ ■ process is radiatively-induced, and the full-one loop calculation was done previously. 

However, for the heavy wino-like or Higgsino-like neutralino, the one-loop cross sec- 
tion violates unitarity, therefore the higher-order corrections may be important. We 
construct a non-relativistic theory for chargino and neutralino two-body states, and 
estimate all-order QED corrections and two-loop corrections by Z and/or W ex- 
change. We find that the critical mass, above that the two-loop contribution is 
larger than one-loop one, is about 8 TeV (O(10) TeV) in the limit where neutralino 
is wino (Higgsino)-like, respectively. Around and above the critical mass, the all- 
order Z and/or W exchange must be included to estimate the cross section. On the 
other hand, the QED corrections depend on the mass difference between the neu- 
tralino and chargino. In the wino-like limit where neutralino is highly degenerate 
with chargino in mass, we find that QED corrections enhance the pair annihilation 
cross section by 1.5-2. 
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I Introduction 



Present observation of cosmological and astrophysical quantities allows a precise 
determination of the mean density of matter (Qm) and baryon density (Qb) m the 
Universe, and the existence of non-baryonic dark matter (DM) is established now 
p. However, the constituent of the DM is still an unresolved problem. The su- 
persymmetric (SUSY) models provides the candidates of the DM, since the lightest 
SUSY particle (LSP) is stable due to the conserved R parity. The LSP may be the 
lightest neutralino in minimal supersymmetric standard model (MSSM). It is a linear 
combination of gauginos (bino and wino) and Higgsinos, which are superpartners of 
gauge and Higgs bosons, respectively. Thermal [2] or non-thermal processes [3] in 
the early Universe may produce the lightest neutralino enough to explain the DM in 
the Universe. 

The detection of exotic cosmic rays is feasible technique to search for the dark 
matter particles, since a pair of the lightest neutralino could annihilate into the SM 
particles with significant cross section |4 . Among those, excess of monochromatic 
gamma ray due to the neutralino annihilation into two photons is a robust signal 
if observed, because the diffused gamma-ray background must have a continuous 
energy spectrum [Sj. 

The neutralino annihilation to two photons is a radiative process. The dominant 
contribution to the cross section comes from the process where a neutralino pair is 
converted into a virtual chargino pair by H^-boson exchange and then the chargino 
pair annihilates into two photons. The full one-loop cross section is calculated in 
Ref. jSj. The surprising fact is that the cross section is suppressed only by the 
W^-boson mass {my/), not by the neutralino mass (m), as ov ~ o^o^jm^^ if the 
lightest-neutralino mass eigenstate is very close to wino or Higgsino. Since other 
pair-annihilation cross sections are proportional to 1/m 2 , the potential of search for 
the monochromatic gamma ray in cases of the wino- or Higgsino-like neutralino DM 
then can be understood as follows; as the neutralino mass increases, the signal rate 
of the monochromatic gamma ray reduces less quickly compared with the other DM 
signatures. 

On the other hand, this cross section should be bounded from above by the 
unitarity limit ov < 47r/(vm 2 ) Thus, for the extreme heavy neutralino, this 
one-loop result must fail jjj, and the higher-order corrections should be included. 
Indeed, the wino- or Higgsino-like neutralino is accompanied with chargino in the 
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same SU(2) multiplet, and their masses are almost degenerate. In these cases, the 
intermediate chargino-pair state in the process is almost on-shell, since the neutralino 
DM is non-relativistic (NR). When m is large, the diagrams are enhanced by a factor 
of a^vnjmw for each W boson exchange, and the higher-order loop diagrams become 
more and more important. Then, the one- loop result is not valid anymore, and we 
need to sum the contributions from ladder diagrams of the weak-boson exchange to 
all orders. 

This failure of the perturbative expansion is similar to the threshold singularity 
in pair creation or annihilation processes in QCD or QED In the quark-pair 
annihilation process at the threshold (NR) region, the higher-order corrections are 
non-negligible since nearly on-shell quark is dominated in the loop integration due 
to the small relative velocity. In other words, the deformation of the quark wave 
function by the QCD potential is not negligible at the threshold region. In order to 
get the reliable cross section, we need to sum ladder diagrams by gluon exchange 
to all orders, or to use the wave function for quarks under the QCD potential. In 
the neutralino annihilation to two photons, we need to include the Yukawa potential 
of the weak boson. When l/(a 2 ^) is larger than the effective range of the Yukawa 
potential 1/mw, the Yukawa potential is point-like in the coordinate space, and does 
not deform the wave functions of neutralino and chargino. However, if m is heavier, 
the wave functions of chargino and neutralino are deviated from plane waves inside 
~ 1/ (a^m) • 111 this case, the perturbative expansion at the threshold region should 
be broken. 

The NR field theory is useful to investigate these higher-order corrections at 
the threshold region in QCD or QED [TJIj. In this technique, we can factorize the 
short-distance physics, such as annihilation or production, from the long-distance 
physics related to the wave function ^T] . The evaluation of the contribution from the 
long-distance physics is possible either by systematical resummation of the ladder 
diagrams in diagrammatic methods or by evaluating the wave function under the 
potential. This formalism is also useful to evaluate the higher-order corrections 
in the neutralino annihilation to two photons. From a viewpoint of the NR field 
theory, the short- distance physics comes from the chargino-pair annihilation to two 
photons, and the long-distance physics from the weak-boson exchanges, in addition 
to the photon exchange between charginos. 

In this paper, the higher-order corrections to the annihilation cross section of the 
neutralino to two photons is studied in cases of the wino- or Higgsino-like neutralino 
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DM. First, we construct the two-body effective action of non-relativistic neutralino 
and chargino pairs. Using this formula, we estimate two-loop corrections by weak- 
boson exchange in a diagrammatic method. The two-loop contribution becomes 
important when the neutralino is as heavy as m^/o^- From the explicit calculation, 
we found that the critical mass, above that the two-loop contribution is larger than 
one-loop one, is about 8 TeV for the wino-like neutralino, and is about O(10) TeV for 
the Higgsino-like one. The resummation of the weak-boson exchange contributions 
to all orders is possible by solving numerically the wave functions of chargino and 
neutralino under the Yukawa potentials of the weak bosons j!2j . We also calculate 
all-order corrections by photon exchange, by using the wave function deformed by 
Coulomb potential. The QED corrections depend on the mass difference between 
neutralino and chargino, 8m. The wino-like neutralino is highly degenerate with the 
chargino, and 5m is typically ~ 0.1 GeV. For this case, we found that the QED 
corrections enhance the cross section by 1.5~2 when the neutralino mass is a few 
TeV. 

This paper is organized as follows. In next section we briefly review masses and 
interactions of neutralino and chargino, especially in the wino- and Higgsino-like 
cases. In section III, we construct the two-body effective action of the NR neutralino 
and chargino pairs and show the strategy for our calculation. Here, in order to show 
the validity for our formalism, we reproduce the previous one-loop result, and show 
that the unitarity bound is satisfied in an extremely heavy neutralino mass limit. In 
section IV, we first include all-order QED corrections to the leading cross section, 
and we evaluate two-loop corrections by Z and W boson exchange. Section V is 
devoted to conclusion and discussion. We give a full effective Lagrangian relevant to 
our study in Appendix A. 

II Non-Relativistic Action of Wino- or Higgsino-like Neutralino and 
Chargino. 

In this section, we review the mass spectrum and the low-energy interaction of 
the SU(2) multiplets containing the lightest neutralino. We are interested in the case 
where the lightest neutralino LSP (xT) is almost degenerate with the lighter chargino 
(x + = xl) m mass. For simplicity, we assume that all the other SUSY particles, 
including the two heavier neutralino states (x° and ^4) an d the heavier chargino 
{xt)i are so heavy that we can ignore their contributions to the DM annihilation. 
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Neutralinos are linear combinations of supersymmetric partners of gauge bosons, 
the bino (B) and the neutral wino (W°), and those of Higgs bosons, the Higgsinos 
(H®, H%). While those four fields have SU(2)xU(l) invariant masses, they are mixed 
with each others by the SU(2)xU(l) symmetry breaking. The lightest neutralino is 
wino-like when M 2 -C \i, Mi and Higgsino-like if \x <C Mi, M 2 . Here Mi and M2 are 
the bino and the wino masses, respectively, and \i is the supersymmetric Higgsino 
mass parameter. When Mi,M 2 ,fi ^> mz, the lightest neutralino mass is 

2 2 

m x\ ^ M 2+ Z Z ° W 2 ^ + /"in 2/3) (1) 
M 2 — /i 

for the wino-like case, and 

777,7(1 + sin 2(3) , „ , 9 , , 9 , 

m *S" " + 20. - W - U) <" - MA - (2) 

for the Higgsino-like case (/i > 0). Here, we show the terms up to 0(m z /msusY)- 

Charginos are linear combinations of the charged wino, W ± , and the charged 
Higgsino (iff, H 2 ). When M 2 ,/i ^> mw, the lighter chargino mass is 



2 

m x+ - M 2 + rJ^ 2 (M 2 + n sin 2/5) (3) 
M 2 — fi 



for the wino-like case, and 



2 

^-^ Ji^n ^ M 2 sin2J) fi) 



for the Higgsino-like case. We show the terms up to 0(m 2 z /uisusy) here, again. 

The mass difference between chargino and the LSP, 8m, is an important pa- 
rameter for the neutralino annihilation cross section into two photons. For the 
wino-like case, their masses are highly degenerate. The tree-level mass difference is 
0{m z jm\ VSY ) in a case m z ,M 2 <C Ml,//, 

4 

Sm tree ~ Ty^s^c^ sin 2 2/3 , (5) 

MifM' 1 

since the 0(m 2 z /msusY) corrections to the masses are SU(2)xU(l) invariant. The 
mass splitting receives the positive contribution from the gauge boson loops due to 
the custodial SU(2) breaking. The radiative mass difference in the wino limit is 

a 2 M 2 ( f m w 2 f ,m z \ 
Sm iad = — (/(— ) - ) - /(0)J , (6) 

where f(a) = Jq 1 dx 2(1 + x) log(x 2 + (1 — x)a 2 ) 13J. This correction is about 0.18GeV 
when M 2 3> m^. 
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For the Higgsino-like LSP, the mass splitting is 0(m 2 z /msusY), 

5m ~ Ym/ w{1 ~ sin2/3) + Ym/ w{1 + sin2/3) • (7) 

The second-lightest neutralino (%°) also degenerates with the LSP and chargino, 
since they are in common SU(2) multiplets. The mass difference between the LSP 
and the second lightest neutralino, 5m N , is again 0(m w j 'msu sy) , 

Sm N = ^f w + ^-s w . (8) 

This is roughly 2 x 5m when tan/? ^> 1. 

Next, we consider interactions among neutralino(s) and chargino. After ignoring 
all SUSY particles except Xi(2) an d X + > the interactions which we have to take into 
account are only the gauge interactions. The leading correction by non-vanishing 5m 
to the neutralino annihilation cross section to two photons is 0(y/ 5m). As will be 
shown in next section, the 0(y/ 5m) correction originates from an infrared behavior 
of the loop integrand which is controlled by 5m. The next-to-leading order correction 
is 0(5m). To this order, the Yukawa interactions of the Higgs bosons or unphysical 
(NG) bosons must be included to keep the gauge invariance. In this paper, we 
calculate the annihilation cross section up to the leading-order correction 0(V5m), 
however, it is straight-forward to include the next-to-leading order calculations. 

In the wino limit, the gauge interactions of inos are 

Ant = -— ffirrWl + h.c.)+e^ Tl^TZ. + eTl^T^. (9) 
sw v * J s w 

The NR Lagrangian can be derived by taking a NR limit of neutralino and chargino 
and integrating out of the gauge fields, and we find 

£ nr } = Vn (idt + ^ Vn + Vc -5m + ^-j rj c + £c (ify -5m- ^-j & 

1 • 



+ f / d 3 y riUx)Zc(v,t) 1— ^ &(y,t)Vc(x) 



\x-y\ 



, . d 3 y Un(x)Vn(M p ^ &(y,t)r ) c(x) + h.c.) . (10) 

zs w J \ \ x ~ y\ ) 

Here, r] C , (,c an d Vn are defined as 

1 + >y 1 — >y 1 + 

Vc = —^Xe- 1 , ic = —^X-e-"* , Vn = A^xV mt . (11) 
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We keep isospin singlet terms only, because only these terms are necessary to calcu- 
late the neutralino annihilation cross section. 

The gauge interactions in the Higgsino limit are 

Ant = (¥irrwi - ^7 M 75x"^ + h.c.) - (\ - c 2 w ) r/r^ 



2sw v ' swcw \2 

+ eF7 M X"^ + 2^^ 15 ^ ■ (12) 
The NR Lagrangian is derived in the same way as the wino-like case, 

45 = ^(^dt+^^m + Vc(^9t-5m+^^r]c 

+ & (id t ~ 6m -^)Zc + & (idt ~ Sm N + ^ £ N 

I — ( 1 ~ 2c M/) 2 e -mz\x—y\ 

+ £ / d 3 y ^(x)e c (y,i) — 4c 7^ ^ 



a 



e -m w \x-y\ 

d 3 y rf N (x)rf N (y,t) — — — £} c {y,t)r) C {x) + h.c. 



8sfy J \ \x-y 

me ( . p -m w \x-y\ , 



d 3 V [tti*)&(M ^ q &{v,t)vc(x)+h.c 



8sw J \ \x-y 



a 



r e -m z \x-y\ \ 



g-m z |rr-j7| 

Here, £ N is defined as £ N = (l—^°)x2^ mt /2, since the sign of the the second-lightest 
neutralino mass is opposite to that of the LSP. 

So far we have derived the NR Lagrangian relevant to the neutralino and chargino 
scattering. Now we would like to include the terms relevant to the neutralino pair 
annihilation into two photons in the framework. The process X1X1 ~^ 77 involves 
external photons whose momentums are of the order of the neutralino mass. Those 
photons cannot be described in the NR Lagrangian. However, it is possible to 
calculate the pair-annihilation cross section in this formalism, by introducing the a 
non-unitary four-Fermi terms in the NR Lagrangian. The optical theorem relates the 
imaginary part of neutralino forward-scattering amplitude T to the pair-annihilation 
cross section a; 

2 Im T = sav . (14) 

Here s and v are square of the center-of-mass energy and the relative velocity, respec- 
tively. The LSP-pair annihilation cross section therefore can be obtained by calculat- 
ing the imaginary part of the forward-scattering amplitude T if proper non-unitary 
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terms are incorporated. This also justifies the factorization of the short- distance 
physics from the the long-distance one, and it is discussed in detail in Ref. [TT] . 

When we explicitly keep the chargino fields in the NR action, the annihilation 
into two photons is described by a chargino four-Fermi operator. The annihilation 
term is given by 

£ ann = d if c (x)£ c (x) tU x )Vc(x) • (15) 

The coupling d is matched to the chargino pair-annihilation cross section to two 
photons using Eq. (|14j). or equivalently va^ n = Im d/2. We find d = irra 2 /m 2 . We 
do not include the annihilation terms for the LSP or the next-lightest neutralino, 
since they are responsible to the processes to W- or Z-boson pairs, not to photons. 



Ill Two-body State Effective Action 

In order to calculate the neutralino annihilation cross section, it is convenient 
to use the effective Lagrangian for the neutralino and the chargino two-body states 
than the NR Lagrangian directly. The two-body state effective action is derived 
by introducing auxiliary fields for ry c ^Ci VnVnj £n£n anc ^ their Hermitian conjugate 
fields. After some calculations, the two-body state effective action is found to be 



S> 



2body 



d A P 

(2^ 



d 3 r 



v a 

25m H 1 h 

m r 




+ E 



8=1 

e 



E - 25m, + 



V 2 



rn 



+ Cc 

2 



r 



4 P)t (r)4 P Vl 



i=l 



-mwT 



-mzT 



+ Cn- 



(16) 



Here, <j)^\r), <j>mif) and 4>^(r) correspond to the chargino pair, the lightest neu- 
tralino pair and the second-lightest neutralino pair, respectively, and 5m± = and 
5mi = 5m n- The capital P is the center of mass energy and momentum of the two- 
body state, and f is the relative coordinate. The internal energy of the neutralino 
pair, E, is defined as E = P° — P 2 / (4m). The coupling constants u)\, u> 2 , (c an d (n 
in Saody are 

0, & = Cv = (17) 



0J2 



Mr 



s 



for the wino-like case, and 

a 



V2a V2a g(l - 2c 2 w ) 2 

4„2 ' ^2 4 2 > 4 2 2 ' ^ 4r 2 2 



for the Higgsino-like case. 

The S-wave states of the neutralino and chargino pairs must give the largest 
contribution to the annihilation cross section, because higher angular-momentum 
modes are suppressed by power(s) of the relative velocity v compared with the S- 
wave state. Thus, the partial wave expansion of £>2body is convenient for us. For this 
purpose, we expand as 

<Vl = E / ^N plm (P)(2p)ji(pr)Y lm (e, 0) • (19) 

l,m 

Here, ji(pr) is the spherical Bessel function, and Yi m (9,4>) is the spherical harmonic 
function. The quantum number p is related to the internal energy E as E = p 2 /m 
under the on-shell condition. In Eq. p9jl . N p i m is the annihilation operator for the 
positive-energy state with quantum numbers (plm). 

As mentioned in Section II, the annihilation cross section of the neutralino pair 
can be calculated from the forward-scattering amplitude by using the optical theo- 
rem, 

m I 2m? " \m 



™(=)=7^**T„(£) (20) 



where the forward-scattering amplitude T PP (E) is 

{2^\P - P') iT pp {E) = ^(0\N p00 (P)Nl 00 (P')\0) (21) 

and the relative velocity v = 2p/m. The amplitude T pp may be obtained by solving 
the equation of motion for the two-body state derived from Eq. (|16|) . as in the 
scattering theory in the quantum mechanics. Or, we may also calculate it by using 
the perturbative theory in a diagrammatic method. 

The dominant contribution to T pp in the each loop order comes from the lad- 
ders of photon and weak-boson exchange, as illustrated in Fig. provided that the 
intermediate neutralinos and charginos are almost on-shell and enhance the corre- 
sponding amplitudes. Similar phenomena are found in QED or QCD as mentioned in 
Section I. In the electron and positron pair annihilation/production at the threshold 
region, the ratio between the lowest-order amplitude and the amplitude with addi- 
tional one-photon exchange between electron and positron is proportional to a/v at 
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Figure 1: Ladder diagrams in the calculation of the forward-scattering amplitude 
{T pp ). The crossing points correspond to the annihilation term for chargino. 

v — > 0. This is well-known as the threshold singularity. In order to evaluate the 
cross section, we need the resummation of the ladder diagrams since a ladder dia- 
gram with n-photon exchange is proportional to (a/v) n . The other diagrams, such 
as the crossed ladder diagrams, are suppressed by additional factors of v. 

The efficient evaluation for the effect of the resummation of the ladder diagrams 
to the electron and positron pair annihilation/production is to use the wave functions 
for electron and positron pair under the QED (Coulomb) potential. When we expand 
the two-body state of the electron and positron pair by the wave functions under 
the Coulomb potential, the Coulomb potential disappears from the two-body state 
action, and the calculation of the annihilation cross section at the threshold region 
is only for a tree-level diagram by the annihilation term, and becomes extremely 
simple. 

In the evaluation for the neutralino annihilation cross section, the resummation 
of weak-boson exchange diagrams, in addition to that of photon exchange, is required 
for the heavy neutralino annihilation cross section. The wave functions of chargino 
and neutralino(s) under both Coulomb and Yukawa potentials can be derived nu- 
merically, not analytically. Then, we also consider the behavior of the cross section 
in a limit of m — > oo. In next section we use the wave functions under the Coulomb 
potential for the chargino pair and the free wave functions for the neutralino pairs 
so that we can incorporate the all-order ladder diagrams of photon exchange, and 
we derive some numerical results. 

In order to demonstrate the validity of the effective action Eq. first, we 

reproduce the previous one-loop result for the neutralino annihilation cross section 
to two photons in our formulation. The leading contribution to Im T pp in the per- 
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turbative calculation is 




2 



Im T pp (0) = 4vra 2 (Anui) 



(2vr) 3 \k\ 2 + m 2 w \k\ 2 /m + 25m 



(22) 



The velocity of the incident neutralinos is almost in the above calculation since 
v/c ~ 10~ 3 for the neutralino DM in our galactic halo. By integrating the r.h.s. of 
Eq. (|22|) and using Eq. flU, we find 



This agrees with the result [Hj which is obtained of the full one-loop calculation in 
the heavy wino or Higgsino limit. 

It is found from Eq. (|22[) that the amplitude is enhanced by a factor of Uim/mw, 
in which m comes from chargino propagators in Eq. (}2*2*|) if 8m can be neglected. 
Other factor 1/m^ comes from the weak-boson propagators and the loop integral. 
The same enhancement factors also appear in the higher-order calculations. When- 
ever one weak-boson exchange is added to the ladder graph, it gives an additional 
factor uim/mw- 



Note that the leading correction due to the mass difference 5m is 0(y 5m) in 

Eq. (J23|) . not 0(5m) as mentioned in the previous section. If the correction was 
0(5m), the integrand in Eq. (jUJ) could be expanded by 5m. However, in this case, 
the integral is infrared- divergent. 

The one-loop cross section is independent of the neutralino mass in small 5m 
limit. On the other hand, the cross section should be bounded by the unitarity 
condition av < 4ir/(vm 2 ). Therefore the higher-order corrections must dominate 
the cross section, so that the correct 1/m 2 behavior would be reproduced. This is 
explained as follow. If the neutralino mass is much larger than the weak-boson mass, 
we can neglect the effect of the weak-boson mass in the equation of motion. This 
can be seen by introducing the dimensionless coordinate x = c^mr in Eq. ()16|). The 
factor e~ mwr can be approximated to 1 in such an extremely large m case. Then, 
it is found that the cross section of the extremely heavy neutralino has the mass 
dependence o ~ 1/m 2 from the dimensional analysis. In fact, we can solve the 
equation of motion and obtain the forward-scattering amplitude analytically in the 
large mass limit, because gauge interactions in this Lagrangian generate a common 
potential 1/r, say the Coulomb force. 

Now we calculate the cross section in this large mass limit. We take the wino 
limit for the LSP as an example, however, the result is similar for the Higgsino-like 




(23) 
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Figure 2: The ladder diagrams of photon exchange relevant to the calculation of the 
forward-scattering amplitude. A is the one-loop amplitude; A ~ aa.2m/mw 

neutralino. It is convenient to change the basis from (0^,0^ ) to ((f) { p,(f) {P) ) m 
the two-body state effective action so that the 1/r potential terms are diagonalized 

as 




Here, A± = 0.5 (a + ( c ± ^J(a + ( c ) 2 + 4u; 2 ). Since A+ > and A_ < 0, 0^ feels 
an attractive force, and 0^ feels a repulsive one. The neutralino annihilation cross 
section is derived from the 0+^ annihilation, while the contribution from 0^ is 
exponentially suppressed due to the repulsive force. Neglecting the 0_ contribution, 
the forward-scattering amplitude is written as follows; 

fan T ^ ~ ^""^ (25) 

Then, the neutralino annihilation cross section to two photons in the large mass limit 
is 

2n 2 ujla 2 \ + 5 1 



m 2 v(A+ + — exp(— 2n\ + /v )) t>m 2 



2.8 x 10" 5 . (26) 



This cross section behaves as av ~ 1/vm 2 , and satisfies the unitarity bound as 
expected. 



IV Higher-order Corrections in Neutralino Annihilation to Two Pho- 
tons 

In the previous section, we showed that the higher-order corrections from the 
weak-boson ladder diagrams must be non-negligible if the neutralino mass m is 
larger than m w /a2- In this section we estimate more carefully the scale where the 
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perturbative approximation is broken down. We also discuss another important 
higher-order corrections from the photon ladder diagrams. 

We start our discussion from the effect of photon exchange. The corresponding 
diagrams in the perturbative calculation of T pp are shown in Fig. El In this figure, the 
orders of the amplitudes are also shown. The lowest-order diagram at the left side of 
the figure gives the amplitude A (A Q ~ aa^mjm-^. The diagrams have an addition 



factor aym/Sm for each photon exchange, and then, the higher-order diagrams of 
photon exchange become important if the mass difference (Sm) is smaller than a 2 m. 



The enhancement factor aym/Sm comes from the massless photon exchange at 
t-channel. For the neutralino DM pair annihilation, the initial neutralinos is almost 
at rest, and then the intermediate charginos just slightly deviate from on-shell states 
by the mass difference between chargino and neutralino, Sm. The situation is similar 
to calculation of the pair production at threshold in QED or QCD. In this case, the 
intermediate particles deviate from on-shell states by 0(v 2 m) where v is the relative 
velocity of the incident particles, and the diagrams have an addition factor a^ s )/v 
for each photon (gluon) exchange. 

If neutralino is highly degenerate with chargino in mass, we need resumma- 
tion of the photon ladder diagrams. Similar to the positron and electron annihi- 
lation/production explained in the previous section, the resummation of the ladder 
diagrams is possible analytically, by expanding the chargino-pair field (0c(r)) i n 
the two-body state effective action Eq. in terms of the eigenf unctions of the 

Schrodinger equation with the Coulomb potential, 

dk 



l,m L' 



/ ^C klm (P)R% l \r) +Y,Cnim(P)R% l \r) 



Y lm {9,</>). (27) 



Here, C^im is the annihilation operator of the continuum state of a chargino pair 
with quantum numbers (klm), and C n i m is one of a chargino-pair bound state with 
(nlm). Rq ' (r) and r!q (r) are eigenfunctions of the radial direction for the states. 
Both Rq and Rc are given by the confluent hypergeometric function, and the S-wave 
parts of these functions are 



TTocmk c _ ikr / am 
exp(— iram/k) \ 2k 



R { ^\r) = 2d- —z~ lkr \F X ( + 1, 2, 2ikr 



R?\rl = 2(^) 3/2 e— /^(-n + 1,2,^) , (28) 

where \F\ is the (Kummer's) confluent hypergeometric function. With this expan- 
sion and Eq. (JTHJ), we derive the S-wave effective Lagrangian of the two-body states. 
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We show details of the Lagrangian in appendix A. 

In the following we calculate the higher-order corrections of gauge-boson ex- 
change. The forward-scattering amplitude T pp is calculated from the S-wave effective 
Lagrangian by treating weak-boson exchange interactions perturbatively. The effect 
of all-order photon-exchange effects are now automatically taken into account in a 
dependent coefficients of the interactions. 

We first calculate the cross section in the leading level where the effect of only 
one VT-boson exchange is included. After some calculation, we find 
'p 2 \ 2txo? 



av 



m I m 2 



iu)\ f dk imAj^c{Pik) I namk 



p J 2ir p 2 — k 2 — 2mbm + i0 + y 1 — exp(— iram/k) 
itoi ^ imA N( j(p,n) ( am\ 3 / 2 



p ^ p 2 + a 2 m 2 /An 2 — 2m5m + iQ~ 
The coefficients A^c and A N q are defined in Appendix A. Here the first term in 
r.h.s. comes from the continuum states of chargino pair, while the second term from 
the discrete states of chargino pair. 

We first investigate contributions of the chargino bound states to the cross sec- 
tion. We can see from Eq. (p?Tj) that their contributions are suppressed by a factor 
a compared with the continuum part. This is because the overlap between the 
wave functions of the neutralino pair and the chargino bound states is suppressed 
by a. If energy of the neutralino pair is on the pole of the chargino bound state 
(p 2 /m ~ 25m — a 2 m/4n 2 ), the cross section is enhanced beyond the suppression 
due to the small overlap of the wave functions. However, it is very unlikely for the 
pair annihilation of the neutralino DM. The average DM velocity in our galactic 
halo is roughly v/c ~ 10~ 3 . Thus, the bound state contribution is negligible as far 
as the neutralino mass is smaller than lOOTeV, since p 2 /m ~ 10~ 6 m and the mass 
difference 8m may be larger than lOOMeV. 

Thus, we only have to take into account the contribution from the continuum 
chargino-pair states. Since the continuum part of the cross section is a smooth 
function around the p ~ 0, the cross section for the neutralino DM annihilation may 
be approximated by av (0). The result is 

av(0) = ^HC(0)| 2 (29) 
m z 



where 

exp 



C(0) = iut 



dk Aim namk 



am 2m\yk 
arctan( — ^ 



2k mlv — k 2 ' 



2?r 1 - e - nam / k k 2 + mfy -k 2 - 2mSm 
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Figure 3: a) The leading XiXi ~^ %1 cross section taking into account all-order QED 
effect, <jv , as a function of the neutralino mass m (solid lines). For comparison, result 
of the one- loop calculation Eq. (J23|) is also shown (dashed lines), b) The contours of 
o"/o"i-ioo P in a (m, 5m) plane. 

In a limit of vanishing a, the av reduces to o"fi_i oop in Eq. 

In Fig. El a), solid lines are av(0) for 5m = O.lGeV, 1 GeV and 10 GeV. The 
left axis corresponds to the wino-like neutralino annihilation cross section, and the 
right axis is for the Higgsino-like one. For comparison, av i_i 00 p in Eq. ()23|) are also 
shown by dashed lines. Note that influences of the QED potential always work to 
enhance the cross section. This is because the Coulomb force acts as an attractive 
force between the chargino pair. In Fig. El b), contours of av(0)/avi-\ oop is shown in 
a (m, 5m) plane. The enhancement is large if the mass difference 5m is small and m 
is large as expected. For the wino-like neutralino, which is highly degenerate with 
chargino as 5m ~ (9(100)MeV, the enhancement by factor of 2 is possible when the 
neutralino mass is (9(l)TeV. 

Next, we discuss the next-to-leading order corrections arising from the diagrams 
of two weak-boson exchange. The detail of the calculation will be given elsewhere 
|12j . and we only show the ratio between the next-to-leading order corrections and 
the leading amplitude, R, in FigJU 

In the wino-like case, the Z exchange between the chargino pair, shown in Fig. 
leads to the next-to-leading order corrections. In Fig. E]a), the ratio R in the wino- 
like case is plotted as a function of m with 5m = 0.1 GeV and 1 GeV. We call the 
value of m where \R\ = 1 as the critical neutralino mass m cr i t . Above and around 
the mass, the higher-order corrections is larger than the leading-order one, therefore 
the perturbation is not reliable any more. For the wino-like case, m CTit ~ lOTeV, 
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Figure 4: The ratio between the next-to-leading order amplitude and the leading- 
order amplitude as a function of the neutralino mass m. a) and b) are for the wino- 
and the Higgsino-like neutralino, respectively. 



Y° 

X l 



Y° 

X i 



X 



X 



w z 

^X ^~ X 



~o ~+ ~ + 

Xi-^-X ^ % 
w z 



X2 



X 



Y° 

X i 



X 



X 



Y°- 

X i 

+ z w 

X l— ^- X2 x 



Wino-like 



Higgsino-like 



Figure 5: Diagrams which contribute to next-to-leading order correction, a) and b) 
are for the wino- and Higgsino-like neutralino DM, respectively. 

and for the Higgsino-like case m crit ~ O(10) TeV. 



V Conclusion and Discussion 

Although the existence of DM is established, we do not yet know the DM na- 
ture. The neutralino DM predicted in MSSM is actively searched for by various 
experiments, trying to solve this myth. The search for gamma-ray signal from the 
center of our galaxy is a feasible way to discover the neutralino DM. Especially the 
monochromatic gamma ray from the neutralino DM pair annihilation is a robust 
signal. 

The reliable estimation of the signal rate is important to either set limit of the 
neutralino DM or interpret the observed signals. The major systematic error of the 
signal rate comes from the distribution of the dark matter in our galaxy. However, 
the pair annihilation cross section cr(XiXi ~^ 2j) nas n °t been also understood 
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completely. The problem arises when Xi is wino- or Higgsino-like. In that case, 
the one-loop cross section crvi_\ oop is approximately constant, and does not scale as 
1/m 2 . The one-loop cross section obviously breaks the unitarity bound in m — > oo 
limit. The enhancement comes from chargino in the loop, which is degenerate with 
the LSP neutralino. 

In this paper we study the higher-order loop corrections to the pair annihilation 
process. We find that the dominant contribution comes from the ladder diagram 
of weak bosons and photon exchange. We study the corrections using NRQED 
technique, which allows us to incorporate ladder-type QED corrections analytically 
to all orders, while including weak boson exchange corrections perturbatively. We 
find the QED corrections enhance the cross section up to a factor of two. We also find 
the critical scale m crit , around and above which the next-to-leading order corrections 
by weak-boson exchange are larger than the leading contribution, is ~ O(10) TeV for 
the Higgsino-like LSP and ~ 10 TeV for the wino-like LSP. The corrections reduce 
the cross section for m < m cr \t- We also find the pair-annihilation cross section 
satisfies the unitarity bound in a limit of m —>■ oo. We note that our formulation 
may be used for any other DM candidate if the DM comes from SU(2) multiplets 
and the DM is degenerate with other components in mass. 

The relic density of the Higgsino- or wino-like LSP with the mass above 1 TeV 
may be consistent to the matter density of the Universe Qm- The monochromatic 
gamma ray in TeV region may be searched for the future atmospheric Cherenkov 
telescopes (ACT) |S] . The signal flux may be above the sensitivity, if the DM density 
distribution is singular at the center of our galaxy. Therefore it is important to 
estimate the neutralino pair-annihilation cross section reliably. For such a heavy 
Higgsino- or wino-like neutralino, one has to sum the weak corrections to all orders 
as we discussed in this paper. The resummation is possible by the solving the NR 
equation of motion of neutralino [T2"] . 
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Appendix A S-wave Lagrangian 
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In this appendix, we list the full S-wave Lagrangian, which is used to calculate 
the leading and the next-to-leading order corrections of the cross section. The S- 
wave Lagrangian is obtained by substituting Eq. (fT^j) and Eq. (|2*7j) to the two-body 
state effective action Eq. (fTEj). and extracting the S-wave parts. We find 



C, 



Cn + Cyj + Cz + C a 



where 



C 



2tt \ m 



cj.c n 



c 



w 



J-y f dpdk . , , . 

i=l n J Z7r 



CtN^ + N^C k 



n p 



r dk 



+ Cv/ 



dp dp' 
2vr 2tt 



c n 



2a 2 



2^ 



namk 



exp(—iram/k) 



\2nJ 



dk' 



amk' 



2tt V 1 - exp(-7ram/A? 7 ) <:7fe ' + V V 2^7/ ^ 



(30) 



Here, we omit the arguments P, I, and m for simplicity, which are the center of 
mass energy and momentum, and orbit-angular momentums, respectively. Coupling 
constants etc, are defined by the Eqs. ()17I18|) . and are annihilation 

operators of the lightest and next-lightest neutralino pairs, respectively, and Cf. and 
C n are for the continuum and discrete chargino pair states. The coefficients A^c 
and A N q, etc, are given by a little complicated functions, 

-iam/2k 

i± i VI (JCIIkK i "'" — " " — '<• ' 

A N c{P,k) = 4/- ; — Im 



am V 1 — exp(— iram/k) 



' m w - i(p + kY 



mw — i(p — k) 
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a i \ 4 /am\ 3 / 2 

A Nc(P> n ) = (tt) Im 

am \ 2n / 



mw — am/2n — ip^ 
vtiyj + am/2n — ip / 



A C c(k,k') 



4 / Tiamk / namk' 



A dd (n, n') 



x 



m z + (k — k') 2 y 1 — exp(— nam/k) y 1 — exp(— iram/k') 

( m z + i{k + k>) Y am/2k f mz + i{k + k') Y aml2k ' 
\m z — i(k — k') J \m z — i(k' — k) J 



( am am 4kk' 
x 2 Fi ( — i + i + 1,2 



2k '2k' ' ' m 2 z + (k - k'Y 
2 a 3 m 3 



x 



4m 2 — a 2 m 2 (l/n — 1/n') 2 (nn') 3 / 2 

1 2m z — am{l/n — l/n')\ n ( 2m z — am{l/n' — l/n)\ n 
\2m z + am{l/n + 1/n') J \2m z + am{l/n + 1/n') y 

4a 2 m 2 /(nn') \ 



x 2 F X 1 - n, 1 - n', 2, 



4m 2 , + a 2 m 2 (l/n — 1/n') 2 



4 / namk /am\ 3 / 2 

c ^ ~ m 2 + (ijfe - am/{2n)) 2 \j 1 - exp(-vram/fc) V 2rT/ 



x 



'm 2 + am/ (2n) + * am ' 2fe ( m z + am/ (2n) + ik\ n 
m z + am/ (2n) — ik ) \m z — am/ (2n) + ik J 



( am 2ikam/n 
x 2 Fi —« + 1,1 -n,2, 



2/c m 2 z — (am / (2n) — ik) 2 ) 

\m 2 + (p — p J' 2 / 
Here, 2 F± is the (Gauss's) hypergeometric function. 
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